Quantum computation, quantum state engineering, and quantum phase transitions 

driven by dissipation 
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We investigate the computational power of creating steady-states of quantum dissipative systems whose evo- 
lution is governed by time-independent and local couplings to a memoryless environment. We show that such 
a model allows for efficient universal quantum computation with the result of the computation encoded in the 
steady state. Due to the purely dissipative nature of the process, this way of doing quantum computation exhibits 
some inherent robustness and defies some of the DiVincenzo criteria for quantum computation. We show that 
there is a natural class of problems that can be solved with such a model — the preparation of ground states of 
frustration free quantum Hamiltonians. This allows for robust and efficient creation of exotic states that exhibit 
features like topological quantum order and the creation of PEPS and it proves the existence of novel dissipative 
phase transitions. In particular the latter can in principle be verified experimentally with present day technology 
such as with optical lattices. 
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The strongest adversary in quantum information science is 
decoherence, which arises due to the coupling of a system 
with its environment. The induced dissipation tends to de- 
stroy and wash out the interesting quantum effects which give 
rise to the power of quantum computation, cryptography, and 
simulation. While this is certainly true for many forms of dis- 
sipation, we show here that dissipation can also have exactly 
the opposite effect: it can be a full-fledged resource for uni- 
versal quantum computation and quantum state engineering 
without any coherent dynamics needed to complement it. 

We consider a quantum system composed of N particles 
(e.g. qubits) interacting with local environments giving rise to 
memoryless and time-independent dissipation processes. We 
will show first how to design the interactions with the envi- 
ronment to implement universal quantum computation. This 
new method, which we refer to as dissipative quantum com- 
putation (DQC), defies some of the standard criteria for quan- 
tum computation since it requires neither state preparation, 
nor unitary dynamics |1|]. However, it is nevertheless as pow- 
erful as standard quantum computation. Then we will show 
that dissipation can be engineered Q] to prepare all ground 
states of frustration free Hamiltonians. Those include matrix 
product states (MPS) ^ and projected entangled-pair states 
(PEPS) jH, like graph states H and Kitaev M and Levin- 
Wenn [7] topological codes. Both DQC and dissipative state 
engineering (DSE) are robust in the sense that, given the dis- 
sipative nature of the process, the system is driven towards 
its steady state independent of the initial state and hence of 
eventual perturbations along the way. 

Apart from novel ways of performing quantum compu- 
tation or state engineering, our results imply that quantum 
phase transitions |8j] may be driven by dissipation alone. That 
is, the physical properties of the steady state of our system 
may change abruptly in the thermodynamic al limit when we 
slightly modify the parameters characterizing the dissipative 



dynamics JjJ. This immediately follows from the existence of 
frustration free quantum Hamiltonians that exhibit quantum 
phase transitions in one [ 1 1 J and two dimensions |l2 ] . 

In this work we will concentrate first on DQC, showing how 
given any quantum circuit one can construct a master equa- 
tion whose steady state is unique, encodes the outcome of 
the circuit, and is reached in polynomical time (with respect 
to the one corresponding to the circuit). Then we will show 
how to construct dissipative processes which drive the system 
to the ground state of any frustration free Hamiltonian. We 
will prove that MPS and certain kinds of PEPS can be effi- 
ciently prepared using this method, with details given in the 
appendix. In this paper we will not consider specific physical 
setups where our ideas can be implemented. Nevertheless, the 
appendix will provide a universal way of engineering the mas- 
ter equations required for DQC and DSE, which can be easily 
adapted to current experiments based on, e.g. atoms in optical 
lattices [13] or trapped ions 1 14]. Thus, we expect that our 
predictions and in particular the existence of quantum phase 
transitions driven by dissipation may be experimentally tested 
in the near future. 

Let us start with DQC by considering TV qubits in a line 
and a quantum circuit specified by a sequence of nearest- 



neighbor qubit operations {Ut\t 



We define \ipt 



U t U t -i...Ui\0)i <g> ...|0)iv, so that \ipr) is the final state af- 
ter the computation. Our goal is to find a master equation 
p = C(p) with Liouvillian in Lindblad form lfl5ll 



C(p) 



k 



L k pL\-\{L\L k ,p} 



(1) 



where the Lk act locally and has a steady state, po : 0) which 
is unique; (ii) that can be reached in a time poly(T); (iii) such 
that ipT can be extracted from it in a time poly(T). As in Feyn- 
man's construction of a quantum simulator 01611 . we consider 
another auxiliary register with states {\t)}J =0 , which will rep- 
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resent the time. We choose the Lindblad operators 

Li = \0)i{l\ ® |0) t <0| (2) 
L t = U t ®\t + l){t\ + U}®\t){t+l\, (3) 

where i = 1, . . . , N and t = 0, . . . , T. It is clear that the L's 
act locally except for the interaction with the extra register, 
which can be made local as well. Furthermore, 



Po 



T + 



^yE I^X^I C4) 



is a steady state, i.e. C(po) = 0. Given such a state, the result 
of the actual quantum computation can be read out with prob- 
ability 1 JT by measuring the time register. In the appendix we 
show that pa is the unique steady state and that the Liouvillian 
has a spectral gap A = ir 2 /(2T + 3) 2 . This means indeed 
that the steady state will be reached in polynomial time in T. 
Note that this gap is independent of N as well as on the actual 
quantum computation which is performed (i.e. independent 
of the U t ). It is also shown that the same gap is retained if the 
clock register is encoded in the unary way proposed by Kitaev 
i30Tl . making the Lindblad operators strictly local. A sketch of 
the proof is as follows: first, we do a similarity transformation 
on C that replaces all gates Ui with the identity gates, showing 
that its spectrum is independent of the actual quantum compu- 
tation. Second, another similarity transformation is done that 
makes C hermitian and block-diagonal. Each block can then 
be diagonalized exactly leading to the claimed gap. 

In some sense, the present formalism can be seen as a ro- 
bust way of doing adiabatic quantum computation 11811 (errors 
do not accumulate and the path does not have to be engineered 
carefully) and implementing quantum random walks 11911 . and 
it might therefore be easier to tackle interesting open ques- 
tions, such as the quantum PCP theorem, in this setting 12011 . 
Also, it seems that the dissipative way of preparing ground 
states is more natural than to use adiabatic time evolution, as 
nature itself prepares them by cooling. 

Let us now turn to DSE and consider again a quantum sys- 
tem with N particles on a lattice in any dimension. We are 
interested in ground states <!/, of Hamiltonians 



(5) 



which are frustration free, meaning that minimizes the en- 
ergy of each H\ individually, and local in the sense that H\ 
acts non-trivially only on a small set A C {1, . . . , N} of sites 
(e.g., nearest neighbors) 12 ill . We can assume the H\ & to 
be projectors and we will denote the orthogonal projectors by 
P\ = 1 — H\. States ^ of the considered form are, e.g., all 



PEPS (including MPS and stabilizer states|22j]). 

We will consider discrete time evolution generated by a 
trace preserving completely positive map (cp-map) instead of 
a master equation. These two approaches are basically equiv- 
alent i23ll as every local cp-map T can be associated to a local 
Liouvillian via C(p) = N\T(p) — p], which leads to the same 



fixed points and a gap that is N times larger|24|]. We choose 
cp-maps of the form 



A 



1 m 



(6) 



where the p\s are probabilities and U\.i, . . . , U\ lW , is a set 
of unitaries acting non-trivially only within region A. They 
effectively rotate part of the high-energy space (with support 
of H\) to the zero-energy space, so that tr[T(p)^] > tr[p^] 
increases. As for Liouvillians (HJ we could similarly take 
L\,i = UiH\, or the ones associated to the cp-map. 

We show now that for every frustration-free Hamiltonian 
the cp-map in Eq.© converges to the ground state space 
if we choose the unitaries U\ i to be completely depolariz- 
ing, i.e., T(p) cx Ea^aP^a + 1a ® tr x [H x p]/\i[l x }. For 
ease of notation we will explain the proof for the case of a 
one-dimensional ring with nearest-neighbor interactions la- 
belled by the first site A = 1,...,N, Assume p is such 
that its expectation value with respect to the projector ^ onto 
the ground state space of H is non-increasing under appli- 
cations of T, i.e., in particular trJp'Jj] = ti[T N (p)^]. Ex- 
pressing this in the Heisenberg picture in which T*{^>) = 
* + EA^Atr A (*)/(d 2 iV) we get 



tr[/>*] > tr[/>*] 
> tr[p*] 



1 



(d 2 N) N 
(d 2 N) N 



■ N N 
. /i=lA=l 



(7) 



where the first inequality comes from discarding (positive) 
terms in the sum and the second one is due to bounding all 
partial traces of H\ from below by the respective smallest 
eigenvalue v. Note that the latter is strictly positive unless 
H has a product state as ground state (in which case the state- 
ment becomes trivial). Hence, we must have ti[pH] = 0, 
i.e., p is a ground state of H. It is easily seen that the same 
argumentation goes through for more general interactions on 
arbitrary lattices. 

The above procedure implies the existence of quantum 
phase transitions driven by dissipation. By changing the pa- 
rameters in the cp-map (or the master equation) one can ob- 
tain that some physical properties of the steady state abruptly 
change, in as much the same way as in Refs. 11 lL 1 1 2fl . 

We have shown that it is possible to engineer dissipa- 
tive processes which prepare ground states of frustration free 
Hamiltonians in steady state. However, in the above proof the 
time for this preparation scales as N , which may be an issue 
for experiment with large number of particles. In the follow- 
ing we give much more efficient method for certain classes of 
frustration free Hamiltonians: commuting Hamiltonians and 
MPS. 

We consider first frustration free Hamiltonians for which 
[H\, Hfj] = and show that the corresponding ground states 
can be prepared in a time that only scales polynomially with 
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the number of particles. The corresponding set of ground 
states contains important families, like stabilizer states (e.g. 
cluster states and topological codes), or certain kinds of PEPS. 
Note that there was no known way of efficient preparation for 
the latter. 

Loosely speaking there are two classes of Hamiltonians of 
this type: (i) Hamiltonians for which all excitations can be 
locally annihilated. In this case the time of convergence scales 
as r = (log TV), (ii) Interactions where excitations have 
to be moved along the lattice before they can annihilate and 
t = 0(N logiV). 

In order to see how the first case can occur consider Eq.© 
and note that it can be interpreted as randomly choosing a re- 
gion A (according to p\ which we may set equal to 1 /N), then 
measuring P\ and applying a correction according the [/'s if 
the outcome was negative. Assume now that when iterating T 
the correction on A does not change the outcome of previous 
measurements on neighboring regions since 



VA^A': [U Kl ,Hy}=0. 



(8) 



In fact, this can always be achieved by regrouping the regions 
into larger ones having an interior I(X) C A on which only 
H\ acts non-trivially |33J and letting the U\ t i solely act on 
/(A). Denote by q the largest probability for obtaining twice 
a negative measurement outcome on the same region A. The 
energy \x[HT M (p)] after AI applications of T decreases then 
as AT(1 - (1 - q)/N) M such that it takes 0({N\ogN)/{l - 
q)) steps to converge to a ground state. The relaxation time 

of the corresponding Liouvillian is thus r = C(log N~). 
Clearly, this is only a reasonable bound if q < 1, a condition 
possibly incompatible with Eq.®. 

Note that for all stabilizer states we can achieve q = 0, since 
there exists always a local unitary (acting on a single qubit) so 
that H\U\H\ = 0. A class of stabilizer states where this is 
compatible with Eq.® are the so-called graph states (jfo. In 
this case A labels (with some abuse of notation) a vertex of a 

graph and H\ = (1 — ai x) U (x ,^£ a ^)/ 2 where aW is a 
Pauli operator acting on site A and £ is the set of edges of the 

graph. Obviously, U\ = af^ does the job. In this special case 
we can get even faster convergence when using the Liouvillian 



(9) 



The corresponding relaxation time can be determined exactly 
by realizing that the spectrum of C equals that of — (H ® 1 + 
l(g)H)/2 so that r = 1 Q. 

Let us now discuss the second type of commuting 
Hamiltonians — those for which Eq.® and q < 1 are incom- 
patible. For this class we can still prove fast convergence 
by making explicit use of the fact that frustration-free ground 
states of commuting Hamiltonians have an efficient PEPS rep- 
resentation. That is, when expanded in computational product 
basis, the coefficients are given by a tensor-network whose 
geometry resembles the lattice structure of the interactions. A 



generic property of PEPS is injectivity [|29|| of local regions 
which is, in fact, a sufficient condition for the state to be the 
unique ground state of its parent Hamiltonian. Consider cases 
of commuting Hamiltonians for which the ground state has 
this property. To specify the cp-map in Eq.© we need to 
sort the regions of interactions Ai, A2, . . . such that the union 
Ajt = Ui=i nas an intersection with Afc+i but does not en- 
tirely cover it, i.e., Ik+i '■= ^k+i\A-k 7^ 0. Such an ordering 
is always achievable by possibly regrouping the interactions 
into slightly larger regions. The reason for this ordering is that 
any U\ k which acts only on 1^, will not alter the energies in all 
regions A; with i < k. That is, we have a weakened version of 
Eq.®. Injectivity of the regions Ik in the PEPS representation 
implies then that there is always a unitary U\ u (a depolarizing 
set would work as well) such that q < 1 [35] and we can thus 
follow the above lines of argumentation. The only difference 
is that due to the weakening of Eq.© the energy does not 
decrease homogeneously, but a low-energy region will grow 
stochastically according to the A^'s which requires extra time 
proportional to the systems size, so that r = 0(N log N). 

There are frustration-free ground states which belong to the 
second class of commuting Hamiltonians but for which injec- 
tivity does not hold (e.g. due to a degenerate ground state 
space). A paradigmatic example is Kitaev's toric code state 
16[] where one has a four-fold degeneracy. Due to lack of in- 
jectivity, we have to prove q < 1 separately which is, however, 
trivial in this case since it is a stabilizer state so that q = 0. 
The action of the cp-map (or respective Liouvillian) can be un- 
derstood as moving all the excitations towards a single point 
where they can mutually annihilate. 

We turn now to another family of ground states of frustra- 
tion free Hamiltonians, namely MPS 13]. Clearly, one pos- 
sible efficient way of preparing them using dissipation is to 
exploit the fact that they can be obtained via a sequential ap- 
plication of quantum gates |25| together with the above DQC 
scheme. In the following we will, however, focus on a differ- 
ent way which does not require an additional time register. 

For the sake of clearness, we will consider translationally 
invariant Hamiltonians, although the analysis can be straight- 
forwardly extended to systems without that symmetry. We 
will specify a cp-map to prepare states of the form 



I*) 



-4, 



Hi ...l N ) 



(10) 



where, the ^4's are D x D matrices. As before, we assume 
the injectivity property which implies that ^ is the unique 
ground state of a nearest neighbor frustration free 'parent' 
Hamiltonian which has a gap. Denoting by p the reduced den- 
sity operator corresponding to particles k and k + 1, and 
Pfc = 1 — Hk will denote the projectors onto its kernel and 
range, respectively. Note that tr(Pk) = D 2 . We take N = 2™ 
for simplicity, but this is clearly not necessary. We construct 
the channel T in several steps. We first define a channel acting 
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on two neighboring particles k, k + 1, as follows 

Hr, e (X) := P k XP k + ^tr(H k X). 

Here, k = 2 r_1 (2c — 1) where r = 1, ...,n and c 

1 2 ,l_r . The action of these maps has a tree structure, 

where the index r indicates the row in the tree, whereas c does 
it for the column. Now we define recursively, 

^r,c • — ~ (^r — 1,2c + <5>r-l,2c+l) + £rR-r,c- (LI) 



Here, r = 2,...,n, c = l,...,2 n ~ r , S hc := K hc , and 
e r+ i = l/M r where M = CW 2 and C> 1 (see appendix). 
Note that S r ,i acts on the first 2 r particles, S r ,2 on the next 2 r , 
and so on. We finally define 



T '■= (1 — en+i)<5 ra ,i + e„+iTZn,2- 



(12) 



In the appendix we show that this map achieves the fixed 
point (up to an exponentially small error in C) in a time 
0(jV lo S2W). The intuition behind the cp-map (Q2) is that 
the channels Si yC , which are the ones that most often applied, 
project the state of every second nearest neighbors onto the 
right subspace. Then S 2 ,c do the same with half of the pairs 
which have not been projected. Then 6>3. c does the same on 
half of the rest, and so on. 

In conclusion, we have investigated the computational 
power of purely dissipative processes, and proven that it is 
equivalent to that of the quantum circuit model of quantum 
computation. We have also shown that dissipative dynam- 
ics can be used to create ground states (like MPS or PEPS) 
of frustration free Hamiltonians of strongly correlated quan- 
tum spin systems. This implies the existence of dissipatively 
driven quantum phase transitions, something which could be 
experimentally tested using atoms or ions. 

Let us stress that we have been concerned here with a 
proof-of-principle demonstration that dissipation provides us 
with an alternative way of carrying out quantum computa- 
tions or state engineering. We believe, however, that much 
more efficient and practical schemes can be developed and 
adapted to specific implementations. We also think that these 
results open up some interesting questions which deserve fur- 
ther investigation. For example, how the use of fault toler- 
ant computations can make our scheme more robust. Or how 
one can design translationally invariant cp-maps that prepare 
MPS more efficiently. Or the importance and generality of 
the set of commuting Hamiltonians, which we believe to be 
intimately connected to the fixed points of the renormaliza- 
tion group transformations on PEPS (as it happens with MPS 
12611). Furthermore, the model of DQC might well lead to 
the construction of new quantum algorithms, as e.g. quantum 
random walks can more easily be formulated within this con- 
text. Finally, other ideas related to this work can be easily ad- 
dressed using the methods introduced, e.g., thermal states of 
commuting Hamiltonians can be engineered using DSE since 



the Metropolis way of sampling over classical spin configura- 
tions can be adopted to the case of commuting operators. Sim- 
ilar techniques could be applied to free fermionic and bosonic 
systems, and, more generally, it should be possible to device 
DSE-schemes converging to the ground- or thermal states of 
frustrated Hamiltonians by combining unitary and dissipative 
dynamics. 
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APPENDIX: EFFICIENCY OF DQC 

In this section we will prove that the Liouvillian defined by 
the Lindblad operators 



U = |0)i(l| ® |0) t <0| 

L t = U t ®\t+l){t\ + U} ®\t){t 



II 



(13) 
(14) 



is gapped. More specifically, it holds that for any initial con- 
dition jo(0), we can show that \\p(t) — p(oo)\\ < e in a time t 
that scales logarithmically in 1/e and quadratically in 1/T. 
The Liouvillian is defined as the (non-symmetric) matrix 

C = ^L a ®L a - g (l^ L l L <* ®I + I®Y1 L « La ) 

a \ a a / 

where a runs over the labels of all Lindblad operators. We 
will bring this Liouvillian into a simpler form by doing two 
eigenvalue-preserving similarity transformations. First, we 
apply the unitary 



W = ^U t U t - 1 ...U 1 ®\t)(t\ 



and observe that the spectrum of the Liouvillian is the same 
as if the quantum circuit would only have consisted of identity 
gates. The spectrum of the Liouvillian is hence independent 
of the actual computation that we want to do (note that this 
was also the case in the context of adiabatic quantum com- 
putation). Without loss of generality, we therefore assume 
Vi : Ut = I. Second, we diagonalize the part acting on the 
logical qubits by the similarity transformation £ = XCX -1 
where 



X = (Jfi®Jf 2 ®..Jfj V )®|00)t(00| + /®(J-|00) t (00|) 
X t - |00)(00| + |00)(11|-|11)(11| + |01)(01| + |10)(10|. 
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Note that the double indices arising in those expressions 
reflect the fact that the Liouvillian acts on the tensor product 
of the physical space with itself. The part acting on the logical 
qubits is now completely diagonal, and C becomes 



C! = \h...ji...)(h...ji...\ ® Li, 

tl ,»2 , . . ,»W ,J1 ,. . .JN=0 

Lij = Z + A i |0) i (0|®7t + A J -7 t ® |0) t (0| 

JV 
k=l 

l = ^ l 4 <g> L t - - ( ]>] L\L t ®i t +i t ®J2 L\L t I 
t V t t / 

L t = \t)(t + l\ + \t+l)(t\ 

The problem of calculating the eigenvalues of L is there- 
fore reduced to the problem of calculating the eigenvalues of 
the T 2 x T 2 matrices for all possible positive integer num- 
bers Xi, Xj. It happens that is block-diagonal, and consists 
entirely of diagonal elements of the form (—2), (—1 — Xi), 
( — 1 — Xj), of 2 x 2 blocks of the form 



-1-Ai 1 W -2 1 
1 -l-Aj'V 2 -2 

and of tridiagonal matrices of the form 



(l + X l + X ] )\0)(Q\-2j:J- 2 1 \t)(t\-\T)(T\ 

+ El=i(\t)(t + M + \t + i)(t\) 

The first classes of diagonal blocks have all eigenvalues 
strictly smaller than —1, and therefore correspond to terms 
that converge extremely fast. The gap is determined by the last 
class of tridiagonal matrices, and that one corresponds exactly 
to the class of matrices that has been extensively studied in the 
context of random walks. The only matrix with eigenvalue 
is the one with A; = Xj = 0, which proves that the fixed point 
is unique. The gap of L o can be calculated exactly, and is 
given by 



1 



r + i 



(T + l) 



The largest eigenvalue for the cases (ij) ^ (00) is obtained 
for Xi + Xj = 1, and is given by 



2 I cos 



2T + 3 



(2T + 3) 5 



It follows that the gap of the Liouvillian is larger than 1/T 2 , 
which we set out to prove. In principle, this does not suffice to 
prove that convergence will be reached at a time in the order of 



1 /T 2 , as the Liouvillian might have exponentially large Jor- 
dan blocks; it can readily be checked that this is not the case 
here. 

Let us next check what happens if we use Kitaev's unary 
encoding J30ll of the time register such as to make all Lindblad 
operators strictly local. We will have to replace L t by 



L t = 7®|l)(l|(g)X. 



I 



where X is acting on the £'th qubit in the unary encoding, and 
add Lindblad terms that force the system to converge into the 
allowed subspace for the time Hamiltonians: 

Lg = J<g>|0)(0|<g>|0)(l|<g>I. 

Furthermore, we have to replace 

Li = (/(g)|0)i(l|®/)®(|0}(0|®/). 

It now happens that the relevant spectrum of the new cor- 
responding Liouvillian is unchanged: all the terms Li, L t , L q 
are such that they map density operators that are in the right 
subspace of the unary encoding into the right subspace, and 
furthermore the only terms that connect the right with the 
wrong subspace are L q which can map wrong states right 
ones. The Liouvillian is therefore block-upper-diagonal, and 
the eigenvalues are completely determined by the diagonal 
blocks. All the eigenvalues in the wrong block turn out to 
be smaller than —1, and therefore the relevant eigenvalue lie 
in the right block. The gap is therefore left unchanged. 



Making use of the recent results of Oliveira and Terhal 113 111 , 
and of Aharonov, Gottesman and Kempe 1531 . it can also 
be shown that the same computational power is retained if 
the Lindblad operators are only acting on nearest neighbour 
qubits in a 2-D square lattice or on nearest neighbours in a 
1-dimensional spin chain of 12- level systems. 

Let us conclude this section by explaining why the gap will 
considerably change if we were to use the master equation 
approach to solve general problems in the class NP. The idea 
would be to put a penalty term on one of the output qubits cor- 
responding to getting the right answer, and relaxing the con- 
straints on some input bits. First note that the construction of 
the Feynman Hamiltonian 11611 is only possible when all quan- 
tum gates are unitary; hence the whole circuit is reversible. 
This implies that the problem only makes sense if some input 
qubits to the quantum circuit are initialized, as otherwise any 
output can be obtained. Because of the fact that both input 
and output qubits must be initialized, we cannot replace the 
actual gates with the identity gate, and as a result the Liou- 
villian does not have a nice block-diagonal structure anymore 
with only blocks of polynomial size, but we get exponentially 
large blocks. Such blocks will typically lead to exponentially 
small gaps. 
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APPENDIX: ENGINEERING DISSIPATION 

Here we show how to engineer the local dissipation which 
gives rise to the master equations (Q]) and cp-maps ©. They 
are composed of local terms, involving few particles (typically 
two), so that we just have to show how to implement those. In 
order to simplify the exposition, we will treat those particles 
as a single one and assume that one has full control over its 
dynamics (e.g. one can apply arbitrary gates). 

Let us start with the cp-maps. It is clear that by apply- 
ing a quantum gate to the particle and a 'fresh' ancilla and 
then tracing the ancilla one can generate any physical action 
(i.e. cp-map) on the system. Furthermore, by repeating the 
same process with short time intervals one can subject the sys- 
tem to an arbitrary time independent master equation. This 
last process may not be efficient. An alternative way works 
as follows. Let us assume that the ancilla is a qubit inter- 
acting with a reservoir such that it fulfills a master equation 
with Liouville operator L a = VTcr_, where er_ = |0)(1|. 
Now, we couple the ancilla to the system with a Hamiltonian 
H = 0(<r_Lt + cr_L). In the limit T 3> il one can adia- 



batically eliminate the level 1 1) of the ancilla 12711 by applying 
second order perturbation theory to the Liouvillian. This is 
done as follows: the unperturbed Liouvillian can be written as 
Co = QDQ^ 1 with Q = Q^ 1 the eigenvectors of Cq\ writing 
the perturbed eigenvectors as Q exp(fi/vTX), we solve the 
equation —XD+DX+QChQ = with Ch the perturbation 
arising from the Hamiltonian part; the effective Liouvillian is 
then given by — fl 2 /TXDX. In this way we obtain a mas- 
ter equation for p describing the system alone, with Liouville 
opeartor n/y/TL. By using several ancillas with Hamiltoni- 
ans H = n(<j-Li + <J_L,\) and following the same procedure 
we obtain the desired master equation. Although we have not 
specified here a physical system, one could use atoms. In that 
case, the ancilla could be an atom itself with |0) and |1) an 
electronic ground and excited level, respectively, so that spon- 
tanous emission gives rise to the dissipation. The coupling 
to the system (other atoms) could be achieved using standard 
ideas used in the implementation of quantum computation us- 
ing those systems 12811 . 



APPENDIX: EFFICIENCY OF DSE FOR MATRIX 
PRODUCT STATES 



We show that the cp-map defined in the text for the creation 
of MPS converges in sub-exponential time. 

Let us denote by 7i r the ground subspace of the Hamilto- 
nian X)I=i Hk corresponding to the zero eigenvalue, and by 
q n the projector onto the orthogonal subspace. The basic idea 
is that by applying the map <S r) i a sufficiently large number of 
times, L r , to any density operator, p, we obtain a state which 
is practically supported on H r - We will show that the error 
fi^ := tr(q r S^ 1 (p) can be made arbitrarily small. In particu- 
lar, for n = log 2 N we just have to choose L = 0(iV Iog2 W). 



In order to simplify the discussion, we will take into ac- 
count the following: S^ +1 1 with M := Le r+ i 3> 1 will be a 
sum of contributions which will typically have the form 



<5 r,l'-V,l ,< "r+l.lO r l O r l l-Cr+1,1 



(15) 



where Li 



appears 

°r-,l'-V,2 



~ (1 — e r +i)/2e r +i and where the channel 7£ rj i 
■> Le r +\ times. We define: L r = l/2e r +i; J r = 

V r+ i = Hr+l,lfr't Pi = /rV r+ \; p[ = T^r+l,\{pl)\ 

Xi = tr(p r+ ipi); Yi = ti(p r +ip'i). We will approximate 

5^ +1)1 ~ Kr+l 1 - Usin § the fact that for k < ^' P kPr = 

PrPk = Pr we have that tr{p r+ i.iS r ^(a)] > ti\p r +i t ia] for 
i = 1,2, and thus Xi > Yi-\. It is easy to show using the 
same thing that 

1 



Yi = Xi + j^tv[ Pr+l tv {QpiQ)]. 



(16) 



Here we have simplified the notation: tr denotes the trace 
with respect to the particles k , k Q + 1, and Q = Qk ,k a +i- 

Now, using that for any projectors P and Q = 1 — P, 
{sfeP ± Q/y/e)p(^feP ± Q/yfe) > it follows that Pl > 

PrPiPr — 3-v/ pr r . Finally, using the properties of MPS 
it is easy to show that tr[p r+ itro(Qp r pzp r Q)]/ D 2 is lower 
bounded by a z(l 
Thus we obtain 



H^ r )(l — X{), where z is a constant. 



X l+ i >Yi>Xi-i^^+z{\ 
Iterating this expression we obtain 



(17) 




-z(l-p^)} 



]ir+l£r+l 



(18) 



After some lengthy algebra, we obtain that if we choose 
e r /e r _|_i = M = CN 2 , and L r = l/e r +i the final error will 
be (3/z) 2 (l — z) c whereas the number of applications of the 
map L\og 2 n = N 2lo ^ N+lo ^ c . Thus, by choosing C suf- 
ficiently large we can always make the error arbitrarly small 
with a subexponential number of applications of the map. 
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